arXiv:hep-ph/9706309v3 21 May 1998 


TIFR/TH/97-12 

UCRHEP-T186 

(revised) 


New interactions in neutrino oscillations with three light flavors 


Ernest Ma[| 

Physics Department 
University of California, Riverside 
CA 92521. 

Probir Roy|^ 

Tata Institute of Fundamental Research 
Mumbai 400 005, India. 


September 1997 

Abstract 


If one assumes solar and LSND neutrino oscillations to explain the corresponding data, 
then the atmospheric neutrino dehcit cannot be accommodated within the Standard Model 
with three light flavors, unless one ignores the data’s zenith-angle dependence. We propose a 
novel solution to this problem by postulating large anomalous diagonal i^T-'quark interactions 
which affect — v-r oscillations traversing the Earth and induce the observed zenith-angle 
dependence. 
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Three flavors of massless lefthanded weakly interacting neutrinos occur in the Standard 
Model. Experimental studies at the Large Electron Positron (LEP) collider have definitively 
established the number of light weak neutrino flavors to be three. However, the masslessness 
of any neutrino is not predicted on fundamental grounds. Beyond the Standard Model, 
theoretical arguments exist showing how neutrinos could acquire tiny Majorana [1] or Dirac 
[2] masses. On the experimental front, there is indirect evidence of small nonvanishing 
neutrino masses from three different kinds of phenomena pertaining to neutrino oscillations. 
(1) The observed depletion [3] of the solar neutrino flux from the prediction of the standard 
solar model in different segments of the solar neutrino energy spectrum, (2) the claimed 
discovery [4] of — z/g and z/^ — z/g oscillations by the Liquid Scintillator Neutrino Detector 
(LSND) experiment and (3) the deficit [5] in the atmospheric neutrino flux, measured on the 
ground in terms of the ratio of ratios R = (z/^ + P^) : (z/g + Pe)expt./ : (^'e + ^e)MC) 
MC standing for the Monte-Carlo expectation - all point to nonzero neutrino masses. The 
canonical best fits to the data from the above three different experimental studies cannot be 
accommodated within the hypothesis of only three light neutrino flavors because of the three 
nonoverlapping ranges of 5m? involved. This and other considerations from astrophysics and 
cosmology have led to the speculation of the existence [6] of a fourth sterile (i.e. electroweak 
singlet) light neutrino as a possible way of reconciling all of the known data. On the other 
hand, if only three neutrino flavors are assumed, two of the three possible neutrino oscillations 
can be explained and the question is to what extent the third can be accommodated. Previous 
studies have chosen the exception to be either the atmospheric data [7] or the solar data [8]. 
A recent general analysis [9] shows that the former hypothesis is in fact favored. 

It seems to us that the results of the solar neutrino and the LSND experiments are quite 
unambiguous, assuming the absence of unknown sources of systematic error. In contrast, the 
detailed conclusions from the atmospheric neutrino experiment seem to depend sensitively 
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on the intricacies of the Monte-Carlo simulations [5] used. It may, therefore, be more prof¬ 
itable to consider three-flavor scenarios which naturally explain the solar neutrino and the 
LSND data and then to explore the observed atmospheric neutrino anomaly within those. 
That leads one naturally to the scheme of Ref. [7]. This scheme, in which the interactions 
are exclusively those of the Standard Model, leads to a universal value of R that can give 
acceptable hts to the sub-GeV and multi-GeV data, integrating over all zenith angles. How¬ 
ever, it disallows any measurable dependence of the data on the zenith angle predicting an 
essentially flat distribution. 

As in Ref. [7], we too start with only three light neutrino flavors and Standard Model 
interactions, but then we extend the latter to include the possibility of anomalous neutrino 
interactions [10]. Specihcally, we allow Vr to have any large nonstandard diagonal four- 
fermion (effectively contact) interaction with quarks. This is motivated by two facts: (1) the 
reported observation [11] of anomalous e’''-quark interactions at HERA, which suggests the 
possibility of anomalous lepton-quark interactions in general; (2) among the three known 
neutrinos, z/^^e are strongly forbidden by experimental constraints to have such large inter¬ 
actions while there exist essentially no restrictions on (In particular, there could be a 
heavy vector boson coupling only to leptons of the third generation but to the light quarks as 
well). As we show below, this will result in the novel possibility of an induced zenith-angle 
dependence for atmospheric neutrino oscillations, which appears to be favored by the data. 

With an anomalous z/T--quark interaction, the survival probability R may vary with the 
zenith angle in atmospheric neutrino oscillations despite a large 5m? chosen at around 0..25 
eV^ to satisfy the LSND data. Since the interaction cross section (in the detector) of neutri¬ 
nos is roughly 3 times that of antineutrinos, and the neutrino flux is somewhat larger than 
the antineutrino flux [12] at higher energies, such a variation is potentially able to explain 
the multi-GeV atmospheric data. The extra Vr interactions inside the sun are offset by a 
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Sm^ much larger than it would be for the canonical matter-enhanced effect, as explained 
below. 

We start with the following approximate mass eigenstates: ui ~ 1^2 ~ Coh'^J. + -So^'r, 

^3 ~ —sot'/i + cqi^t, where cq = cos 6 * 0 , sq = sin 6 * 0 , and 9q is not small. We choose mi ~ 0, 
m 2 ~ 10“^ eV, and m 3 ~ 0.5 eV. We then allow ui to mix with z /3 with a small angle 6' and 
the new ui to mix with z /2 with a small angle 6. The exact mass eigenstates are then 

Z/i = cc'ue + Cs'{-SoUf, + CoUr) - s{CoU^ + SoUr), ( 1 ) 

z /2 = c(CoZ/^ + SoZ/r) + ScVe Ss'(-SoZ/^ + CqZ/^), (2) 

= c'(-SoZ/^ + CQUr) - s'Ve- (3) 

Let us now turn one-by-one to the three sets of neutrino oscillation data. 

LSND; ~ oscillations, as probed in this experiment, are controlled by ~ 0.25 
eV^. Any signihcantly higher value chosen for 5m|^ will be in contradiction with restrictions 
imposed by the search for z/^ disappearance in the CDHS experiment [13] for large angles 
(which will be needed later in explaining the atmospheric neutrino effect). On the other 
hand, for such a value of 5m?, the LSND data [4] imply a mixing angle y with sin 2% ~ 0.19. 
These numbers are just about compatible with the constraints of the Bugey experiment [14]. 
Comparing with (3), we hnd 2sos'c' ~ 0.19. 

Solar neutrino data; The canonical solution for solar neutrino oscillations takes vi ~ z/g 
and z /2 ~ a linear combination of z/^ and Vr, with m 2 > mi and some mixing between z/i 
and z/ 2 . In its passage through the sun, z/g gets an extra induced mass because of its forward 
scattering with the electrons. The matching of this mass with 5m\^ produces the well-known 
MSW effect [15]. Here we have new extra diagonal z/T--quark interactions. Consequently, a 
larger 5m\^ is needed to cancel against the induced Vr mass, which should be negative in 
this case [16]. 
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Let us examine in detail the effect of z/r-quark interactions on the passage of electron 
neutrinos through the sun. We write this new interaction as [10] 

/Inew = -V2iyrat,^rL{Gl^vQl^q + <^rrA97''75g) (4) 

for all quarks q. Note that the r neutrino and the r antineutrino are known from r-decay 
properties to be lefthanded and righthanded respectively. Therefore, all possible four-fermion 
interactions involving them and quarks can be brought into the form (4) by Fierz transfor¬ 
mations. Only the vector coupling contributes to the potential relevant to forward 

scattering while its contributions for the neutrino and antineutrino cases are equal in mag¬ 
nitude but opposite in sign. We dehne = Gl^y/Gp as in Ref. [10]. Note that we do not 
require any flavor-changing interactions which would have been necessary to obtain oscilla¬ 
tions if Sm'^ = 0. 

Solar neutrino oscillations occur between and z/q, = cqz/^ + SqUt- with angle 9. Since 
ml 3> ml, Up = —Squ^ + cqu^ effectively decouples. Hence the relevant evolution equation 
can be written, after rotating away a common phase, as [10] 

4,E,-(‘'A c, ( ° V"']. (5) 

" dt \UaJ \mlsm29 2 ml cos 29 + 4:y/2GFE,^{sle'^Ng — N^) J \i^a / 

In (5), Nqdq = and Ne^u,d is the number of (electrons, u-quarks, d-quarks) per 

unit solar volume. The coefficient Sq = sin^ 9q of the d^Ng term in the second diagonal matrix 
element originates from the 3 x 3 —2 x 2 flavor matrix reduction. In order to have a large 
Cg and yet satisfy the resonance condition for solar-neutrino flavor conversion, we see that 
m 2 should be larger than its canonical value of 2.45 x 10“^ eV, and dg should be negative. [If 
dg comes from R-parity violating squark exchange, then it must be positive; but if it comes 
from vector exchange, then it may be of either sign.] 

We now assume as a crude approximation that Ng ~ 4iVe in the sun. It then follows from 
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(5) that the effective mixing angle for z/g — Va oscillations in solar matter is given by 


tan26*£, = 


sin 29 


cos 29 + 2 v^G_pm 2 — l)Ne 

and the MSW resonance condition [15] is 


( 6 ) 


m 2 + 1 ) = cos 29{2V2GFNeE^)-\ (7) 

In the canonical MSW solution, the left-hand side is (6 x 10“®eV^)“^. It can thus be matched 
with the requirement 


~ -3.92 = -4.17(m2/10-W) + 0.25. (8) 

Our seemingly arbitrary choice of Srri^i ~ 10“^ eV^ is now seen as a reasonable value so 
that e'g can be large enough to be relevant for the following discussion on the atmospheric 
neutrino data. Note that the same range of 9, i.e. near 0.04, works here as well as in the 
standard MSW solution for solar neutrino oscillations. 

Atmospheric neutrino data: The depletion in the flux of muon neutrinos and antineu¬ 
trinos, produced in the earth’s upper atmosphere, is caused by flavor oscillations 

which occur between the physical states z /2 and z /3 with 6 m ‘^2 ~ 0-25 eV^ and angle 9o- The 
oscillation wavelength “in vacuo” is A ~ AnEy/ 5 ml 2 ~ 10(i7;y/GeV) km, Eij being the neu¬ 
trino energy. Hence, for E^ < 10 GeV, several oscillations occur in the earth’s atmosphere. 
Gonsequently, one obtains the classical survival probability Pq = Cq -|- Sq = 1 — |sin^26*o. 
The choice of Sq — 0.47 yields P = Pq ~ 0.66 for the ratio of ratios, in reasonable agreement 
with the data for downward going neutrinos and antineutrinos. This implies s' ~ 0.21. 

For neutrinos and antineutrinos coming downward, the density of the atmosphere is 
negligible for the new diagonal z/T--quark interactions to be of any importance. For upward 
moving ones, the density of the earth turns out to be in the right range for them to make a 
difference. Since z/g gets effectively decoupled from the z/^ — Vr oscillation problem, we now 
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have 


AiE, 


dt V 


m| sin 26Q 


m- 


sin 26*0 2ml cos 29o + Ay/^GpE^e'^Nq J \u. 


( 9 ) 


For the earth we estimate an average iVg ~ 9 x 10'*^ Thus if one chooses to dehne 

the parameter X = e'qE^/{10 GeV), the effective mixing angles in terrestrial matter (a) 6^ 
between and Ur and (b) between z/^ and 9^ are respectively given by 

sin 26*0 , sin 26*0 


tan 20^ COS20O + O.O91X 


, tan 20^ COS20O-O.O91X’ 


( 10 ) 


with cos 26*0 = T0.091X as the resonance conditons. 


For sub-GeV neutrinos, the X term is insignificant, but for multi-GeV neutrinos it may 
become large enough for the resonance condition to be satished. Assuming adiabaticity, the 
neutrino and antineutrino flavor survival probabilities are described well by the formulae [17] 


PK- 

= ^(1 + cos 26*0 cos 26*,^), 

( 11 ) 

P{i>,- 

^9fj) = ^(1 + cos 26*0 cos 20 ,^), 

( 12 ) 


where cos 26*,^ and cos 20,^ are computed from ( 11 ). Although the conditions for adiabaticity 
may not be satished, our purpose is to try to hnd the maximum effect for a given magnitude 
of the anomalous interaction which is of course unknown. Any nonadiabaticity would only 
tend to diminish this effect. 


Owing to the opposite signs of the media contributions to neutrino and antineutrino 
oscillations, matter effects will get somewhat diluted. However, there are two important 
factors to be considered. First, the initial z/^ hux is larger than the z/^ hux for multi-GeV 
neutrinos. In the upper atmosphere, ^^{9^) is produced together with p+(p“) from 7 r’''( 7 r“) 
decay. The subsequent decay of fj,~^{ij,~) to 9^{i/^) will equalize the total number of z/^ and 
9^, but there is an energy dependence and given that the jjG /ratio is larger than one [ 12 ], 
we allow a factor of r = ratio of the z/^ to 9^ hux in our following discussion. Second, in the 
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detection of atmospheric neutrinos, there is no experimental measurement of the charge of 
the resulting leptons. Specihcally, /x-like events include both /i“ and jj,'^, but they are not 
separated. Now ~ Sup for an isoscalar target, hence the measured probability is weighted: 


3rP + P 
3r + 1 


(13) 


In the absence of media effects, P = P even with oscillations (assuming CP conservation), 
hence P^ = P = P. In the presence of media effects, P 7 ^ P, so the above expression for 
Pm should be used. 


Suppose we now make the same choice of Sq — 0.47 which leads to a suppression proba¬ 
bility Pq = 0.66 for the downward travelling neutrinos. Then, for an optimistically large and 
negative value of X = —15, we get P = 0.31 and P = 0.76, hence Pm is lowered to 0.39 if 
r = 1.5 in Eq. (14) or 0.42 if r = 1.0. This is a potential explanation of the observation of 
a smaller R for atmospheric neutrinos going through the earth (zenith angle > 7 r/ 2 ) where 
R = Pm than for those coming down through only the atmosphere (zenith angle < 7 r/ 2 ) 
where P = Pq. 


There have been studies [18] of P-parity violating squark interactions, scalar and vector 
leptoquarks, as well as contact interactions of neutrinos. Whereas these may be related to 
our proposed effective interaction, they are all restricted to be small. We consider instead 
a vector boson B which couples to u'y^^u -|- with coupling qb- This 

interaction would result in a negative X as desired. Using the result of a previous model [19] 
where B couples to baryon number, and allowing for the fact that here B couples to u and 
d but not s, c, and h, we hnd as = fi's/dvr < 0.057. We now require Qb/^'b — 15Gir, hence 
= 64 GeV is allowed. The deviation from e — fi — r universality in Z decay is then given 
by 

(i"i + i"2) = ^(0.5735)(1.32) = 0.0069, (14) 

^T^\9L + 9RJr 27r 
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where Fi and F 2 are well-known functions [19] of the ratio m^/M^ which is about 0.5 for 
= 64 GeV. The above is exactly two standard deviations from the experimental data 
[20], taking into account the kinematical correction due to nir in Z ^ r+r“. The deviation 
in the total invisible width of Z is 0.012 versus the standard-model value of 3, which is again 
two standard deviations from the data [20]. This is then a possible explicit model for our 
scenario. 

Our scenario will have the following consequences for the forthcoming experiments. Solar 
neutrino experiments will conhrm the MSW solution, but the interpretation of is subject 
to the ambiguity that it could be Sm^/(1 — d^ge^) instead. However, the anomalous z/^-quark 
interactions will be observable at SNO, thereby resolving this ambiguity. Both z/g and 
G conversion experiments will measure a at around 0.25 eV^, but sin^ ~ 0.036 
while sin^ ~ 0.69. The former is outside the region being probed by reactor experiments 
(z/g disappearance) such as Chooz and Palo Verde; the latter is outside that being probed 
by short-baseline accelerator experiments (z/,- appearance) such as CHORUS and NOMAD. 
On the other hand, both regions are covered by all the proposed long-baseline experiments 
(either through z/^ disappearance, or z/g and i/j- appearances) such as MINOS, K2K (KEK- 
PS/Super-Kamiokande) and CERN-SPS/ICARUS. 

More immediately, the new data from Super-Kamiokande, Soudan 2, and MACRO on 
events through the earth will be sensitive to the anomalous z/,--quark interactions. 
There should be an energy dependence as well as a zenith-angle dependence. In particular, 
the zenith-angle dependence should be absent or much smaller for sub-GeV data. For a zenith 
angle near zero, our proposal is easily distinguishable from the ~ 10“^ eV^ oscillation 
interpretation because we have R = Pq whereas the latter would require i? ~ 1, owing to 
the short distance between production and detection in that case. To test our hypothesis 
further, the detection and acceptance efficiencies of neutrinos versus antineutrinos have to be 
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understood in more detail. Better yet, the capability of these experiments for distinguishing 
neutrinos from antineutrinos should be explored [21], 
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